Combined approach of density functional theory and quantum Monte Carlo method 
to electron correlation in dilute magnetic semiconductors 



OO 
O 

o 

(N 

o 

Q 



O 

CZ3 



Jun-ichiro Ohc 1 , Yoshihiro Tomoda 1 , Nejat Bulut 1 ' 2 , Ryotaro Arita 3 , Kazuma Nakamura 3 , 

and Sadamichi Maekawa 1,2 
1 Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 
2 CREST, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan 
3 Department of Applied Physics, University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan 

(Dated: December 2, 2008) 

We present a realistic study for electronic and magnetic properties in dilute magnetic semicon- 
ductor (Ga,Mn)As. A multi-orbital Haldane- Anderson model parameterized by density-functional 
calculations is presented and solved with the Hirsch-Fye quantum Monte Carlo algorithm. Results 
well reproduce experimental results in the dilute limit. When the chemical potential is located be- 
tween the top of the valence band and an impurity bound state, a long-range ferromagnetic correla- 
tions between the impurities, mediated by antiferromagnetic impurity-host couplings, are drastically 
developed. We observe an anisotropic character in local density of states at the impurity-bound- 
state energy, which is consistent with the STM measurements. The presented combined approach 
thus offers a firm starting point for realistic calculations of the various family of dilute magnetic 
semiconductors. 

PACS numbers: 75.50.Pp, 75.30.Hx, 75.40.Mg, 71.55.-i 



C 

O 

o 



> 
o 

CO 

o 

OO 

o 



X 



The discovery of ferromagnetism in dilute mag- 
netic semiconductor (DMS) materials, represented by a 
(Ga,Mn)As system, has created much activity in the field 
of spintronics [2, 0> H, El- I n the low-doping regime 
(<C 1%), (Ga,Mn)As is insulating with a clear experi- 
mental evidence for the presence of an Mn-induced im- 
purity band located at 110 meV above the valence band 
[5| . The position of the impurity level and the impurity- 
induced carriers are key quantities in generatin g th e fer- 
romagnetism. Various model studies [& 0, H, SHH have 
clarified that this ferromagnetic correlations between the 
impurities are based on the polarization mechanism of 
the impurity-induced carriers. 

The challenge for predicting new functional properties 
of the DMS materials and their realistic designs clearly 
requires many tasks such as nonempirical electronic-band 
structure of the host compound, reliable estimations for 
impurity-host hybridizations, and rigorous treatments 
for local electron correlations at the impurity site. A 
combined approach of density functional theory (DFT) 
[ll| and quantum Monte Carlo (QMC) technique [T3| is 
clearly useful for this purpose. The hybridization or mix- 
ing parameters between the host and impurity orbitals 
is a crucial parameter characterizing material differences 
and therefore this parameter should be estimated on the 
basis of ab initio density- functional calculations. On the 
other hand, density-functional calculations often fail to 
describe local electron correlations, so this problem must 
be treated with more accurate solvers such as the QMC 
technique. Critical comparisons with the experimental 
data are needed for a check whether this approach is re- 
ally reliable in describing the details of the "low-energy" 
electronic and magnetic properties and/or whether it has 
an extended applicability for many DMS materials. 

In this letter, we present a comprehensive study for the 
generation mechanism of the ferromagnetism of a dilute 



(Ga,Mn]As system. A multi-orbital Haldane-Anderson 
model [13| is employed to describe electronic/magnetic 
properties of this system. First principles band-structure 
calculations are performed to evaluate the mixing matrix 
elements in this model and the resulting model is solved 
with the help of the QMC technique to examine the inter- 
impurity and impurity-host magnetic correlations. The 
hybridization between host and impurity electrons gen- 
erates the impurity bound state jl4j |. When the chemi- 
cal potential is located between the impurity bound state 
and the valence-band top, a notable long-range ferromag- 
netic correlation between the impurities and its dramatic 
enhancement as the temperature decreases are observed. 
We will show that this ferromagnetic correlation is me- 
diated by antiferromagnetic impurity-host couplings. In 
addition, our calculated host local density of states is 
highly anisotropic around the impurity sites, which is 



-■-T=730 K 




" T=360 K 




-•-T=180 K 










/// 







-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 
H (eV) 

FIG. 1: (Color online) Square of the impurity magnetic mo- 
ment ((M z ) 2 } vs the chemical potential fi at various temper- 
atures in the single-impurity case. These results are for the 
single-impurity Haldane-Anderson model. The vertical line 
((j, = 0) denotes the top of the valence band. 
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FIG. 2: (Color online) Inter-impurity magnetic correlation FIG - 3; ( Color online ) (M'M2) vs. |R|/o at T = 360 K for 
function (MfMf ) vs. the impurity separation [R|/o at vari- various values of fj,. 
ous temperatures for fi fixed at 0.05 eV. 



also consistent with the experimental STM image [15j . 
These quantitative agreements with the experiments in- 
dicates that the presented combination idea is a useful 
approach for aiming at realistic calculations and search- 
ing new functional high-T c DMS materials. 

A multi-orbital Haldane- Anderson model is given as 
follows [3: 

U= ]T (e kQ - mKIc£ q + £ {E di - n)d$<% 

k,a,i,£,cr z,£ 

Here, c£„ (c£ a ) creates (annihilates) a host electron with 
wavevector k and spin a in the valence or conduction 
bands denoted by a, and (d^) is a creation (anni- 
hilation) operator for a localized electron at impurity d 
orbital £ located at transition- metal site i. The first term 
in Eq. ([1]) represents a kinetic energy ek Q of a host ka 
electron, and the second term describes a bare onsite 
energy of the impurity and fi is a chemical poten- 
tial. The third term specifies impurity-host hybridization 
Vka.;£ and the last term represents the onsite repulsion 
U at the impurity site. In the present study, we keep 
only diagonal terms in the onsite interaction by neglect- 
ing the off diagonal terms in the interaction (i.e., Hund's 
rule coupling). 

When a Mn atom is substituted in place of a Ga site, 
the five 3d orbitals of the Mn ion are split into the three- 
fold degenerate t 2g orbitals and the two-fold degenerate 
e g orbitals by the tetrahedral crystal field. In the present 
Mn 2+ case with a 3d 5 configuration, the e g orbitals are 
fully occupied and therefore are inactive for the magnetic 
properties. In contrast, the three-fold ti g orbitals being 
a partially-filled state are active for the low-energy mag- 
netic excitations. 

The onsite repulsion was set to £7=4.0 eV following a 
photoemission spectrum [l6[ and this value is very close 
to theoretical values 4-5 eV for metallic Mn [TtJ ■ We note 
that both U and E^ can shift the energy of the impurity 
bound state. Experimentally, it is widely agreed that an 



impurity band exists in the dilute limit of (Ga,Mn)As 
at 110 meV above the valence band ||. To realize this 
situation, we found that E^ should be set to —2.0 eV. 

Estimations for the hybridization parameters Vka,i£ in 
Eq. (fi} need careful treatments. Empirical estimations 
with the Slater-Koster table fail to describe the disper- 
sion of the conduction bands, thus leading to a serious 
failure in the hybridization-parameter estimation. We es- 
timate the hybridization parameter Vka.jf in Eq. |TJ from 
ab initio density- functional calculations [181 ] . The matrix 
element is given by 

Vka,i£ = (lpka\'Ho\tPi£) , ( 2 ) 

where Ho is the one- body part of Eq. ©, \ fa a ) = c£>), 
and \if.^) = g&'|0). \t/fk a ) can be expressed with a linear 
combination of Wannier functions as 

i^> = ^EE a ^(fe)e 4kR kRM)> (3) 

* fi R 

where |u>r m ) is the Wannier orbital at the lattice R 
and the label /x specifies the type of the orbitals. In 
the present GaAs case, the \i specifies the four types of 
sp 3 orbitals centered at a Ga or As site. Expansion co- 
efficients {a afi (k)} are obtained from solving the equa- 
tion Tlcrylipka) = ekalV'ka)- By substituting Eq. © into 
Eq. @, we obtain the form of 

V ka<ii = -^EE<M( fe ) e_tkR (^l^o|^>. (4) 

fi R 

For an actual evaluation of the matrix element 
(wB. l j,\l-lo\ l Pi£ )j we utilize a supercell calculation: We first 
arrange a supercell containing 3x3x3 fee primitive cells, 
where a Ga atom located at the origin in the supercell 
is replaced with an Mn atom. We then make maxi- 
mally localized Wannier functions [l!| for this system 
[2fl |. The width of an energy window was set in or- 
der to include all the valence-bonding and conduction- 
anti-bonding states. The total number of the resulting 
Wannier functions are 217 in which there are 212 GaAs 
Wannier orbitals {|uVi^)} with fj, and n specifying four 
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types of the sp 3 orbitals and sites of Ga or As, respec- 
tively, and five Mn-3<i Wannier orbitals {|u>^)} sitting 
on the origin site i. Then, we calculate the matrix el- 



ement \W n u,\H 



3x3x3| 
KS 



"it, 



with TC 



3x3x3 
KS 



being the Kohn- 



Sham Hamiltonian for the supercell. This matrix element 
is used instead of (wr^I^o!^) in Eq. (j4|) for R ^ i. We 
note that the error due to this replacement is scaled by -tt 
with M being the total number of the atoms in the super- 
cell. For R = i, (wR l i\'H.o\<Pi£) is expected to be negligibly 
small, since wave-function symmetry of |u>r m ) and 
are different. We found that the hybridization parame- 
ters are rather short-ranged; these are nearly zero except 
for the nearest and next-nearest sites; for example, the 
parameters between the Mn-J 2s orbitals and the nearest 
As-sp 3 orbitals are 0.16-0.74 eV, while the values between 
the t2 g orbitals and the next-nearest Ga-sp 3 orbitals are 
0.06-0.14 eV. Notice that, if the empirical Slater-Koster 
estimation is employed, the latter next-nearest values are 
zero. 

We perform the quantum Monte Carlo simulation for 
a one- or two- impurity Haldane- Anderson models, in 
which we use the Hirsch-Fye QMC algorithm [l2| for ob- 
taining the electronic correlations in the one- and two- 
impurity Haldane- Anderson models. In particular, we 
study with QMC the magnetic correlations and the lo- 
cal density of states around the impurity site. We de- 
fine the magnetization operator at the £ orbital of the 
i'th impurity site as Mf^ = — n^i, and at site r 
of the host as m z (r) = U|(r) — ^j.(r). Here, and 
n CT (r) are the spin-dependent number operators at the 
impurity and host sites, respectively. For the single- 
impurity case, we present data on the average value of 
the square of the local moment at the impurity site, 
((AF) 2 ) = jEe=i(( M ^ 2 ); wh ere £ sums over the t 2g 
orbitals, while for the two-impurity case we show data on 
the magnetic correlation function defined by (M-fMJ) = 
| ^^ =1 (M^M| ? ) for impurities located at sites Ri and 
R2. In the single-impurity case, we will also present 
data on the impurity-host magnetic correlation function, 
(M z m z (r)} = |EjLi( M (X TO2 ( r ))> where r is the host 
site and the impurity is located at the origin. 

Figure [T] shows the square of the impurity magnetic 
moment ((M z ) 2 ) as a function of the chemical poten- 
tial fj, at various temperatures T for the single-impurity 
Anderson model. With decreasing T, we observe that a 
step-like discontinuity develops centered at y, « 0.1 eV, 
which indicates that the impurity bound state is located 
at this energy [10l| . Within the Hartree-Fock picture 
S fiot , the impurity- host hybridization creates both the 
impurity bound state (IBS) and the host split-off state at 
luibs- For (i < cjibS) the hybridization with the valence 
band reduces the impurity magnetic moment. When /1 is 
increased to above wibSj the polarization of the split-off 
state cancels the polarization of the valence band, which 
in turn decreases the screening of the impurity moment. 
This is the reason for why ((M z ) 2 ) increases at wtbs- 

Figure O shows the impurity-impurity magnetic corre- 
lation function (M-fMf) vs |R|/a at various T for the 
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FIG. 4: (Color online) Impurity-host magnetic correlation 
function (M z m z (r)) vs |r|/o for [s, = 0.05 eV and various 
temperatures in the single-impurity case. (Inset) Local den- 
sity of states of the host electrons in the (011) plane. Here, 
the impurity is located at the next layer, and the circles in- 
dicate the position of the Ga and As sites. The size of this 
figure is 5 x 5.6 nm 2 , and these results are for T = 360 K and 
H = 0.05 eV. 



two-impurity Anderson model. Here, |R| = |Ri — R2I is 
the impurity reparation and a is the GaAs lattice con- 
stant. In addition, [i is fixed at 0.05 eV so that it is 
located between the top of the valence band and the im- 
purity bound state. In Fig. [21 we observe that the ferro- 
magnetic correlations become stronger and extended in 
real space, as T decreases down to 180 K, which is in the 
physically relevant temperature regime for (Ga,Mn)As 
Q. Figure [3] shows (MfM|) vs \R\/a at 360 K for 
/j = -0.15, 0.05, 0.15 eV. For fi = -0.15 eV, the chemi- 
cal potential is located in the valence band. In this case, 
the system is metallic and we observe Ruderman-Kittel- 
Kasuya-Yosida (RKKY) type oscillations in the magnetic 
correlations. The ferromagnetic coupling disappears for 
/i = 0.15 eV, where the impurity bound state becomes 
occupied. These results indicate that the chemical po- 
tential should be located between the impurity bound 
state and the valence band for obtaining extended ferro- 
magnetic coupling between the impurity moments. This 
case corresponds to the insulating state of (Ga,Mn)As 
found in the dilute limit Q. In addition, the strong n, 
R and T dependence of (MfMf) observed in Figs. [3] and 
U makes it difficult to determine phenomenologically the 
magnetic coupling in DMS materials, instead requiring 
approaches which treat the band structure and correla- 
tion effect on an equal footing and exactly such as the 
DFT+QMC method used here. 

In order to determine the origin of the ferromag- 
netic correlation between the impurities, we present in 
Fig. [U the impurity-host magnetic correlation function 
(M z m z (r)} vs \r\/a at various T for the single-impurity 
Haldane- Anderson model. These results are for /1 = 
0.05 eV in which case the inter-impurity ferromagnetic 
correlation is strong. We observe that the impurity- 
host magnetic coupling is anti-ferromagnetic and it en- 
hances with decreasing temperature. When the impu- 
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rity bound state is occupied (/i = 0.15 eV), we find that 
the impurity-host anti-ferromagnetic magnetic coupling 
becomes much weaker (not shown here) and the inter- 
impurity ferromagnetic correlations collapse. Hence, the 
anti-ferromagnetic impurity-host coupling generates the 
indirect ferromagnetic coupling between the impurities. 

Finally, in the inset of Fig. [U we show local den- 
sity of states around the Mn site defined by N(r) = 

£<n(r) - n(oo)) 



, where n(r) is the electron 

M=0.05eV 

number operator at site r. The highly anisotropic distri- 
bution of local density of states around Mn is consistent 
with the STM measurements [l5|, [22j , which again em- 
phasizes the necessity of inclu ding both the band struc- 
ture and correlation effect [HI, |24| . 

In summary, we have investigated the magnetic prop- 
erties of the dilute magnetic semiconductors (Ga,Mn)As 
in the dilute limit within a combined approach of the 
density functional theory and the quantum Monte Carlo 
method. For this purpose, we have mapped the elec- 
tronic state of the system to the multi-orbital Haldane- 
Anderson model by using the density functional theory 
with maximally-localized Wannier functions. For calcu- 
lating the magnetic properties, we have used the Hirsch- 
Fye quantum Monte Carlo method, which does not ex- 



hibit a fermion sign problem. Within this approach, we 
have correctly reproduced the experimental value for the 
energy of the impurity band of (Ga,Mn)As in the dilute 
limit. The numerical results emphasize the role of the im- 
purity bound state in producing the ferromagnetic corre- 
lations. We have found that, when the chemical potential 
is between the top of the valence band and the impu- 
rity bound state, the inter-impurity ferromagnetic corre- 
lation originating from the anti-ferromagnetic impurity- 
host coupling is strong and extended in real space. This 
case corresponds to the insulating state found in the low- 
doping regime of (Ga,Mn)As. In addition, we have ob- 
served a highly anisotropic distribution of local density of 
states around Mn which is consistent with the STM mea- 
surements. The agreements found with the experimental 
data on the dilute limit of (Ga,Mn)As suggest that the 
proposed approach is a reliable method for investigating 
and designing other DMS materials. 
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